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Abstract 

In this work is discussed possibility and actuality of Lagrangian ap- 
proach to quantum computations. Finite-dimensional Hilbert spaces used 
in this area provide some challenge for such consideration. The model dis- 
cussed here can be considered as an analogue of Weyl quantization of field 
theory via path integral in L. D. Faddeev's approach. Weyl quantization 
is possible to use also in finite-dimensional case, and some formulas may 
be simply rewritten with change of integrals to finite sums. On the other 
hand, there are specific difficulties relevant to finite case. 



1 Introduction 

One initial reason to consider possibility of Lagrangian formalism in quantum 
computations was an idea to use analogues of a minimal action principle in 
some quantum optimization algorithms (c/ Ref. P). 

In quantum physics the minimal action principle has standard correspon- 
dence with classical one: any possible trajectory V of some particle has con- 
tribution exp(-| J^C{t)dt), where C is Lagrange function. Such an expression 
has an oscillatory behavior, and only a trajectory near extremum of action 
A = C{t)dt matters, because other paths compensate each other due to in- 
terference (c/ Ref. 121). 

Does it possible to use such principle in the quantum information science 
for finding an extremum of some function using a quantum mechanical system 
with appropriate Lagrangian? The idea encounters specific obstacles because 
Hilbert spaces used in quantum computations are finite-dimensional. On the 
other hand some analogue of Lagrangian theory may be really built and the 
present paper describes some basic properties of such models. In Section El are 
briefly revisited ideas and formulas used in continuous case, and necessary for 
a further revision with finite-dimensional Hilbert spaces discussed in Section 13 
A relevance to the theory of quantum computations is discussed in Section 01 
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2 Continuous case revisited 



2.1 Weyl quantization 

Let us recall briefly the idea of Weyl quantization [2j with most attention to 
topics necessary for further applications to finite-dimensional Hilbert spaces. 

In continuous case Weyl quantization uses [31^ a function f{p, q) with two 
real arguments p, q and with Fourier co-image /(a, P) described by expression 



fip, ^^=JJ exp(iap + i(3q)f{a, (3) da d(3. (1) 
Such a function is associated with an operator / defined as 

/ = / / exp(mp -f i/3q)f{a, f}) da d(i. (2) 



Here p and q are momentum and coordinate operators. In exponential form used 
in Eq. Q it is defined a Weyl system, i.e., two families of operators: 

lJ[a) = exp(mp), V{I3) = exp(i/3g) (3) 

satisfying relation ^ (in system of units with h = 1) 

U{a)V{l3) = eyi^{ia(3)V{l3)U{a), (4) 

equivalent to usual Heisenberg relation for coordinate and momentum 

[v, q]^pq-qp = -i- (5) 

Formally Eq. |0J follows from Eq. ^ due to a CampbcU-Hausdorff formula 
exp(a + b)~ exp(a) exp(6) exp(— - [a, b]) exp i?2, 

where R2 is a term with commutators of higher order, i?2(p, q) = 0. Using this 
formula it is possible also to rewrite Eq. jSJ as 

f = JJ f{a,(3)exp{ia(3/2)U{a)V{(3)dadp. (6) 

Weyl quantization via Eq. Q and Eq. ^ produces a method of construction 
of some operator / for any function f{p, q) with two variables W : / — *■ /. It 
is enough to use the inverse Fourier transform in accordance with Eq. to 
create / and use it for construction of the operator /. 

There is also direct formula for /. Using an expression for kernel of Eq. 12Jl 
in coordinate representation 0] 

(qY-P+^^''^\q") = cxp(z^-t^/3) Siq" - q' + a), (7) 
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it is possible to write elements of the operator / 0] 



{<l'\fW') = ^J f[P^ ' ^dp. (8) 

It is simple to invert Fourier transformation on p used in Eq. ^ and using 
notation q = {q" + q')/2, x = q" ~ q' to write 

/(p, q)= J {q- x/2\f\q + x/2) exp{~ipx) dx. (9) 

The formula Eq. lO show, how to calculate W^^. 

It should be mentioned, that for a particular case, then / is the statistical 
operator ("density matrix") p, Eq. coincides with an expression for Wigner 
function. The property is known also as Wigner- Weyl isomorphism JS^ . Despite 
of such coincidence in mathematical expressions, Weyl quantization has rather 
different area of applications. Wigner function was suggested for presentation 
of mixed states, but Weyl quantization often used for description of evolution of 
arbitrary operator in Heisenberg picture without actual necessity to work with 
states and density matrices. 

The Hamiltonian function here H{p, q) — W^^H displays a proper corre- 
spondence principle with the classical physics, i.e., the evolution law for opera- 
tors in the Heisenberg picture applied to coordinate and momentum 

| = | = (10) 

produces canonical Hamiltonian equations for classical coordinate, momentum 
and function H{p, q) [2| 

dq ^ dH{p, q) dp ^ dH{p, q) 

dt dp ' dt dq ' ^ ' 

In continuous case such property ensures correspondence principle, because any 
operator may be presented formally via series with p and q. In the discrete case 
there are no good analogues of such operators and so other ideas should be used 
instead. 



2.2 Lagrangian approach and path integral 

The functional (path) integral in Weyl approach may be introduced by quite 
straightforward procedure @j. In Schrodinger picture \ip') = S'j^/') and it is 
possible to write elements of evolution operator S ("quantum gate") for small 
period of time At as 

{q"\SiAt)\q') = {q"\e-'"^'\q') = W'l^ " tHAt\q'). (12) 
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Applying Eq. © to Eq. (O 



27r 7 

It is possible to divide finite interval of time on N small periods and use expres- 
sion S = n exp(i_ff(t)At). For N oo and At — > such products of Eq. ifT^ 
corresponds to expression for functional (continual, path) integral along all paths 
in coordinate and momentum space 

{q"\SW) = l^cxp{iJ{pq-Hip,q))dt)VpVq. (14) 

The function A = J{pq — H{p,q)^dt here is limit At — > of finite sums in 
Eq. 11311 and coincides with classical action. More details about finite sums 
used in the passage to the continuous limit may be found also in description of 
Eq. ^ below. 



3 Finite-dimensional case 
3.1 Weyl pair 



An advantage of Weyl relations Eq. Q for purposes of present paper is possi- 
bility to write an analogue of such operators in finite, 7i-dimensional case 3 . It 
is enough to use Weyl pair of n x n matrices IJ, V: 

2ni 

Ujk = (5j,fe+l mod n, Vjk = exp(— — j)5jfe (15) 

with property 

UV = exp(27n;/ri)V"C/. (16) 

Let us consider finite analogues of other expressions. Let M(a, h) is function 
of two integer arguments a, 6 = 0, . . . , n— 1. It is possible to use discrete Fourier 
transform for both arguments and write analogue of Eq. Q 



1 /o ■ \ 

M(p, g) = - V exp —{ap + bq) M{a, h). (17) 

Tt ^ ^ \ T? / 



Ji-1 

n 

a,b=0 

The function M(p, q) is defined for any real p, but integer values are enough 
to find M{a,b) via inverse Fourier transform, and so value of M{p,q) for any 
real values. 

Analogue of operator / defined by Eq. © is n x n matrix M could be written 

as 

1 l\ / \ 
M = - V exp ( — a6 I JJ'^y^Af (a, 6), (18) 

a, 6=0 ^ ^ 
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but it produces some problems due to using of modular arithmetics, e.g. unlike 
with continuous case the operator Eq. H18|l M is not Hermitian for real function 
M{p,q) Eq. ((T7|l . Anyway, for simpler illustration of suggested approach, it is 
enough to consider asymmetric construction, i.e., 



1 

M' = - ^ U^V''M{a,b), (19) 

a,h=0 

For such operator true discrete asymmetric analogue of expression Eq. ||HJ| 

n-l 



1 " 

(j|Af'|fc) = -J2 e^^^^~'''^M{p, k). (20) 
It can be checked directly using Eq. H17|) and Eq. H19|) : 

n—l 1 

i y e'^P('=-^)M(p, fc) = i y e'^(f('=-^)+'^P+'"=)Af(a, 6) 

p—O p,a,b—0 

- n—l n — l n—l 

= i 5^e^'"=M(a,6)^e^^('=-^+'') ^ ^ <5,_,^,e^'"= M(a, 6) = 

a, 6=0 p=0 a,b=Q^ Z ''^ ' 

^ p ' {U'^v^, 

ndj-k,a 

Inverse transformation for Eq. H2U|I also may be simply found 

n-l 

M{p,q) = Y,e~"-^P=Ml^^j, p,qeZ. (21) 

It should be mentioned also, that in equations above notation for modular arith- 
metic is often omitted for simplicity, e.g., q + j used instead of q + j mod n. 

3.2 Lagrangian approach to discrete models 

Using Eq. H2U|) it is possible to write analogue of Lagrangian function Eq. 

n-l 



{j\S^\k) - l^e^^'(^-'=)(l-zi/(p,j)AO 

^ n — l 

^ L^(.^PU-k)-tH{p,j)At_ ^22) 



p=0 

Here {j\S\k) = Skj are simply indexes of matrix. For interval of time divided 
on N segments we have 

n-l 

{qo\S\qN) = X! {qolSilqi) ■ ■ ■ {qN-i\SN\qN)- (23) 

gi,...,gjv-i=0 
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Using Eq. it is possible to write 



71 — 1 

{qo\S^\q^)^— ^MKp) (24) 

qi,...,qN-l=0 
PO;---,Pn -1=0 

where discrete analogue of action 

^9,p ^ Y ( — Pk{<lk+i - Qk) - H{pk, qk)At) 

N-1 . 

= Yi^P'^^~H{pk,qk))At (25) 

k=0 

is calculated along all n^~^ possible paths qt between two fixed points and 
with different momentum pk for each separate segment of such broken line. 
Approximate exprcsion Eq. H24|l converges^ to value 5'^^^^^^ in limit N oo, 
At —> 0. For continuous limit n oo the Eq. (|25(l for asymmetric ordering may 
look even more traditional [S]. 



3.3 Precise expression with "effective Lagrangian" 

For further analysis of discrete models it is possible to introduce some "effective 
Lagrangian" to have precise value of Sjf. for any division of time interval, includ- 
ing = 1. Really, errors in expressions above was related with two consequent 
approximations like exp(ie) = 1 + ie. Let us instead of application of (discrete) 
Weyl quantisation to Hamiltonian apply it directly to operator S*' in Eq. H22I) . 

1 ""^ 

(Jl^'lfc) = -5]e^P(^-'=)5(p,j) 

^ n—1 

= L'^e^pU-k)+inS{p,j) ^26) 

^ p=0 

with straightforward generalization to sum like Eq. (|24|) . precise for any N ^ 1. 
It is only necessary to change H to H'^^ = « In S'(Ai)/At, it is an "effective 
Hamiltonian," H'^^ ^ H, At ^ 0. Here all paths with terms S{p, q) — should 
be simply omitted. 



3.4 Discrete Fourier transform 

For some cases the precise expression Eq. (|26|) may be quite simple. Let us 
consider for example discrete Fourier transform 

^Strictly speaking, such limit for path integrals is usually not well-defined. 
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Using inverse transformation Eq. 121|l it is possible to find 



F{p,q) - 




Here "effective Hamiltonian" was rewritten in symmetric form {q^ +p^)/2 
due to term Spq. It is convenient also, because for quantum computation with 
continuous variables the Hamiltonian of harmonic oscillator really may be used 
for realization of Fourier transform [Jj . 

4 Discussion 

In general, quantum computation is theory about efficient solving of compu- 
tationally hard problems using quantum systems and processes. Difficulty of 
computational problems relevant with Lagrangian approach and calculation of 
sums over paths has two reasons. First one is usual exponential growth of 
Hilbert space dimension with respect to number of quantum systems, i.e., cal- 
culations with n qubits requires operations with 2" x 2" matrices (multiplica- 
tions, exponents, etc.). The second one is well-known difficulty with definition 
and calculation of functional (path) integrals, it may be described as unspecified 
growth of number of terms in sums like Eq. ^ due to using expressions with 
limits At 0, N ^ oo. 

It was already mentioned some mathematical resemblance with theory of 
Wigner functions in continuous case. For discrete models some methods from 
theory of Wigner function like doubling of lattice, using reflection operator 
and Galois fields |H1 El do not have simple extensions for theory considered 
here, but maybe it provides promising challenge. Anyway, it is quite likely, 
that due to some analogy between theory of Weyl quantization and Wigner 
function, it is possible to combine such techniques for particular tasks related 
with statistical description of ensembles in theory of quantum computations and 
communications . 

In considered theory Fourier transform is related with simplest quadratic 
Hamiltonian similar with harmonic oscillator. In usual, continuous theory only 
quadratic Hamiltonians produce possibility of more or less rigor calculation of 
path integral. In discrete theory used here all expressions may be calculated 
for any functional dependence, but simple representation of Fourier transfor- 
mation may be convenient, because it is a fundamental tool for many quantum 
algorithms ^01 including Shor's factoring one 

It is also possible to compare such approach with yet another propositions 
for quantum algorithms, despite they may look different from traditional one. 
For example it was suggested |12) to use many-slits interference with different 
paths for resolution of NP-complete problems on quantum computer. 

It should be mentioned also idea of quantum optimization algorithms dis- 
cussed in introduction. It works, if due to oscillatory behaviour and appropriate 
Lagrangian, all paths except optimal are vanishing. It is usual for classical limit. 
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On the other hand, the method may be useful for calculation of sum itself, 
even if it is not related directly with any optimization task. For example, ap- 
plication of sum over paths for computing of polynomial equation over finite 
fields was discussed recently in |T3]- It was also noted there, that many basic 
papers about quantum computing complexity use some variants of sum over 
path approach (see ^2] and references therein). 

Some difference between the just mentioned models and the method con- 
sidered in this paper similar with distinction between earlier version with sum 
other spatial paths and more recent version with paths in the phase (p, q) space 

The models discussed in |13| exploit straightforward representation of matrix 
multiplication via path summation Eq. (|23() with only one kind of indexes (qk), 
but in present paper is also used more general expressions with sum other paths 
in whole discrete phase space {pk,qk) Eq. f^ . Only for some particular cases, 
like quadratic Hamiltonians or "effective Hamiltonian" with d symbol Eq. (|28|l . 
the summation on pk may be cancelled. 

The detailed representation in phase space is essential, because one purpose 
of this paper was to find relations between physical models and abstract quan- 
tum computations, e.g., it becomes clearer from equations considered above, 
how Lagrangian term pAq — HAt appears in expressions for sum over paths, 
how classical limit may be obtained for quantum system due to growth of Hilbert 
space dimension, etc. 

There are also other tasks relevant to given approach. For example, La- 
grangian formalism is standard method for description of interacting quantum 
fields and it may be useful for more adequate models of quantum gates and 
computation. 
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